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We study the dynamics of quantum coherence and quantum correlations in two semiconductor double-dot 
molecules separated by a distance and indirectly coupled via a transmission line resonator. Dominant dissipation 
processes are considered. The numerical results show the sudden death of entanglement and the robustness of 
quantum discord to sudden death. Furthermore, the results indicate the dephasing processes in our model can 
lead in the revival and decay of coherence and discord with the absence of entanglement for certain initial states. 
By observing the dynamics of coherence versus discord for different initial states, we find that the similarities 
and differences of coherence and discord are not only related to the dependance of discord on optimizing the 
measurement set, but more importantly to the coherences in individual qubits which are captured by the adopted 
coherence measure. 
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I. INTRODUCTION 

Quantum entanglement, as one of the most intriguing hall- 
marks in the quantum world, manifests nonclassical correla- 
tion between quantum systems which cannot be found in any 
classical system UJJ] . Entanglement is a key ingredient for 
quantum communication, quantum cryptography, and quan- 
tum computing [2]. Besides the entanglement, other non- 
classical correlations have been demonstrated 0431- Ollivier 
and Zurek introduced the so called quantum discord (QD) 
to describe the difference between two quantum expres- 
sions of the classical mutual information. Nonclassical corre- 
lation may still exist with the absence of entanglement, which 
implies the quantum discord as a more proper measure of 
the quantumness of correlations. When quantum discord is 
vanishing, the local accessibility of the classical information 
without perturbing the composite system is shown. Recently, 
quantum discord has received great attention on various top- 
ics, e.g., its potential to be the resource for deterministic quan- 
tum computation with one quantum bit |6, 7], the immunity to 
the "entanglement sudden death" (ESD) phenomenon 18 , |9 1, 
the sufficient condition for completely positive m aps ifloL [Till , 
the QD close to the quantum phase transition Ill2l - fl4f1 , the be- 
havior of QD in the Grover search [15], the monogamic rela- 
tion between the entanglement of formation (EoF) and the QD 
iflrjl . and the relation between entanglement irreversibility and 
the OP lfl7ll. 

On the other hand, among a number of physical systems 
UH], spins in semiconductor double-dot molecules (DDMs) 
are promising candidates for qubits to implement quantum 
information processing (QIP) in a solid-state system. Ap- 
proaches based on DDMs combines spin and charge manip- 
ulation and thus take advantage of the stability lfl9ll . the scal- 
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ability [20], and the efficiency for readout and coherent ma- 
nipulation of charge states [21]. Several schemes adopting the 
architecture with DDMs have been proposed to show poten- 
tial for low-noise coherent electrical control [22], to generate 
cluster states 112311 and to implement quantum computing B24I1 . 
Moreover, Fanchini et al. M25I1 investigate the dynamics of 
EoF and QD between two DDMs, considering the effects of 
dissipation. In their work, each DDM has one excess electron 
and the two DDMs are coupled through direct Coulomb inter- 
action, which is fine for demonstrating the quantum correla- 
tions between two solid-state qubits. However, if the direct in- 
teraction between qubits is strong, it may be problematic when 
a third adjacent qubit is introduced. The direct interactions be- 
tween the third qubit and the former two qubits may disturb 
the operation of the two-qubit gates and the control of non- 
adjacent qubits by indirect interaction may avoid this. This 
motivates us to investigating the dynamics of quantum cor- 
relation between nonadjacent qubits, and verify whether the 
considered quantum correlations, QD and EoF, are resistant 
to dissipations in a scalable solid system. For this purpose, 
we adopt a quite different architecture: each DDM contains 
two excess electrons and the two DDMs are indirectly cou- 
pled through a transmission line resonator (TLR) I126l - l28ll via 
a capacitor. Dissipation processes, e.g. the photon leakage, 
the pure dephasing and energy relaxation of qubits are taken 
into account for realistic. 

Furthermore, we are interested in the relation of quantum 
correlations and quantum coherence (QC). Quantum correla- 
tion and coherence both arise from quantum superposition. 
If there exists coherence among subsystems of a compound 
quantum system, the nonlocal coherence may induce the non- 
classical correlation between the corresponding subsystems, 
besides the local coherence of each corresponding subsystem 
l29tl . In this case, quantum correlation is a kind of nonlocal 
coherence. We investigate the coherence dynamics from the 
point of off-diagonal elements of the density matrix describ- 
ing the composite system. The effect of pure dephasing to the 
dynamics of QC and quantum correlations are presented and 
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demonstrated. The similar and different behaviors of QC and 
QD are shown and discussed for varies initial states. 

The paper is organized as follows. In Sec. II, we briefly 
review the basic concepts of quantum discord. In Sec. Ill, 
we describe the adopted architecture in detail and give the 
Hamiltonian of the system. The dominant dissipation pro- 
cesses, the typical system parameters, and the master equa- 
tion are demonstrated. In Sec. IV, we present and discuss the 
numerical results for the effect of the pure dephasing and the 
behaviors of quantum coherence and quantum discord in two 
coupled DDMs. Finally a summary and some prospects are 
provided in section V. 



II. QUANTUM DISCORD 

In classical information theory, the correlation between two 
systems J{ and S can be measured by the mutual information, 
: S) = 7f(^)+W(£)-7fG?l,£), where "HQ denotes the 
Shannon entropy ^Hip) = - Yijk Pjk ^og 2 Pjk 021] with the prob- 
ability distribution p. The first extension of the expression to 
the quantum case is to directly replace the Shannon entropy 
by the von Neumann entropy Sip^'s) - -Tr(pj7!s log 2 p^s), 
where pjvs is the density operator describing the bipartite 
quantum system. The quantum mutual information between 
systems J[ and S is thus defined as 

?(pm) = S (p$d + S (ps) - S (pwed, (1) 

where p^ = Trgp^g, p<g = Tr^p^g is the reduced density ma- 
trix obtained by taking the trace over all states of the system 
S or M, respectively. 

Through the Bayes rule, the classical mutual information 
can be equivalently rewritten as J($l : S) = ( H(Jl)- ( H(Jl\S), 
where ( H{J[\'B) is the conditional entropy which quantifies the 
ignorance about the state of J[ on condition that the state of 
S is determined. When generalizing to the quantum case, the 
conditional entropy depends on what measurement is locally 
performed on T> and what outcome is obtained. If one fo- 
cuses on projective measurements described by a complete set 
of one-dimensional projectors {Ma }, the remaining state of 
based on the outcome k is p^ — (I ® Mk)pyis(I ® M k )/p k , 
where I is the identity operator performed on J[ and p k = 
Tr(I ® Mjt)p^s(I ® M k ) is the probability for obtaining k. The 
quantum conditional entropy is thus defined as the the residual 
entropy, 5 (p^|{Ma}) = Y*kPkS{p:n\k)- The second extension 
of the classical mutual information to the quantum case can 
be written as the difference of entropy of J{ before and after 
the measurement on S: 

J(pms){M k ) = S (p^ - ^ p k S ip nk ). (2) 

k 

All the nonclassical correlations between J[ and £ are re- 
moved by the local orthogonal projective measurement on 
S while contained in both the initial and residual entropies, 
therefore ^J(pff\ r B){M k ) only captures the classical correlation 
between subsystems. S(pji\{M k }) depends on the choice of 
{M k }. To maximize SF(fiym){Mu} one should optimize the set 



{M k } to seek the minimum of S (p^\{M k }). Therefore the clas- 
sical correlation is defined as 14J 

C(pj\<b) = S (pa) - min (Mt | ^ PkS (pn\k)- (3) 

k 

The discrepancy between the two quantum extensions of 
classical mutual information defines the quantum discord JH]: 

QHpas) = Hpsm) ~ C{p mi} )- (4) 

For a bipartite quantum system, which is just the case in this 
paper, Hamieh et al. Bill have demonstrated that the pro- 
jective measurement is the positive operator valued measure 
which can maximize the classical correlations. This property 
is useful in our numerical calculation process. 

III. DESCRIPTION OF THE MODEL 

The architecture we adopt consists of two identical DDMs 
capacitively coupled to a TLR, as shown in Fig. 1 . Each DDM 
is formed by a GaAs/AlGaAs heterostructure with a layer 
of two-dimensional electron gas below. Several lithograph- 
ically defined metallic gates are used to control the double- 
well potential for two dots and interdot tunneling [21]. With 
a modest external magnetic field B z = lOOmT along axis 
z, Zeeman effect results in energy gaps between the spin- 
aligned states |(1, l)T+> = ITT), 1(1, l)7-> = III) and the spin- 
anti-aligned states |(1, 1)T > = (IU> + lit))/ V2, |(1, 1)5) = 
(IT4) _ lit))/ V2. The notation (ni,n r ) labels n/ electrons lo- 
calized in the left dot and n r electrons in the right dot. Follow 
the approximation of Ref. JUIHH, each DDM can be de- 
duced to a two-level system in the basis {|(1, 1)5) , |(0, 2)S)}, 
described by the Hamiltonian 

H d = -A |(0, 2)S ) <(0, 2)S | + T (|(1 , l)S ) <(0, 2)S | + h.c.) (5) 

where H = 1, and A and T are the relative potential and tun- 
neling [19] between |(0,2)S) and |(1, 1)5), respectively. The 
energy gap between the eigenstates is a>d = V4r 2 + A 2 . Here 
we choose the optimal point A = for maximum coupling 
JH, and use the resulting |1> = (|(0,2)S) + |(1, 1)S»/ V2and 
|0) = (|(0,2)5) - |(1, 1)5))/ V2 as qubit states. 

We consider the TLR with length L, capacitance per unit 
length Co, and characteristic impedance Zq. The DDM is 
on resonance with the fundamental frequency of TLR, idq = 
7t/LC[)Zq, therefore we can neglect all the higher energy 
modes of the TLR j20h . The Hamiltonian of the TLR reads 

H r = ojoa^a, (6) 

where a (a^) is the annihilation (creation) operator of the res- 
onator. 

As shown in Fig. 1, the two DDMs are individually coupled 
to the TLR via the identical capacitor C c . The charging energy 
provides the DDM-TLR interaction, H mt - C c V r Vd, where V,- 
and Vd are the voltages of the TLR and DDM, respectively. 
The Coulomb interaction between the two DDMs can be ne- 
glected due to the long distance (~ 4p.m) between them with 
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FIG. 1. (Color online) Schematic diagram for the adopted architec- 
ture. The TLR has a length of L, the DDMs has a relative potential A 
between two singlet states. The two DDMs are individually coupled 
two the TLR via the identical capacitor C c . The TLR is connected to 
the input/output ports via a capacitor C e for writing/reading signals. 



respect to the small size of DDM (~ 400nm). After the quan- 
tization of the volta ges o f the TLR [26], V r = Vw/(LCo)(a f + 
a) and the DDM Hi, V d = 2e/C lol |(0, 2)5 ><(0, 2)5 1 + 
e/C to ,\(l, 1)5)((1, 1)5 1, and in the rotating wave approxima- 
tion, the full Hamiltonian takes the usual Jaynes-Cummings 
form 

H = i ^ <o di o- zi + w a 1 a +8^ ( fl0 ~; + hc ^ W 

i=l,2 (=1,2 

where cr z is the pauli z operator, cr'" = |1) (0|, cr = |0) (1|, and 
the coupling coefficient 

= eC c I gjp (g) 
2C /0/ \ LCo 

I 



Here C tot is the total capacitance of the DDM. As mentioned 
above, we consider the two DDMs identical and resonant with 
the fundamental frequency of TLR, u^x = todi - <^o, thus the 
Hamiltonian in the interaction picture reads 

V = gYj(a°-] +h-c). (9) 

i=l,2 

Now we demonstrate the dominant dissipation processes 
considering the interaction with environment in the Born- 
Markov approximation. First, the photon leakage of the res- 
onator which can be described by the damping rate k — u>o/Q, 
where Q is the quality factor of the TLR. Though the inter- 
nal loss of the resonator is negligible (Q m t up to 10 6 ) i26ll27ll . 
with the external magnetic field B : the quality factor will be 
reduced to Q ~ 10 4 13211 . The TLR frequency chosen as 
a)o/2n = lOGhz leads to the damping rate A-/27T = u>o/27tQ = 
1 Mhz. The second is the spin relaxation induced by the qubits 
coupling to the phonon bath. The induced spin relaxation time 
is obatined as 7\ ~ 1 ps using the spin-boson model at the 
optimal point II22I1 . The third is the spin dephasing results 
from the low frequency fluctuations of the electrostatic bias 
and the hyperfine interaction with the nuclear spins. The for- 
mer dephasing time is estimated as T2 ~ tOdT^ ~ 10 ns ll22ll . 
with the observed bare dephasing time Tb ~ 1 ns j33ll . The 
latter time-ensemble-averaged dephasing time is measured as 
~ 10 ns and by spin-echo technique the dephasing time 
may be prolonged beyond 1 /is 1I2II1 . Taking into account the 
above processes, the evolution of the state of the three party 
system (the two DDMs plus the TLR) can be described by the 
master equation lf28ll 



p = -i[V,p] + V -^(<r zi pcr Z i -P)+Y\ ■^■(2o- 1 pcr! - cr) cr,p - per] cr,) + ^(2apa 1 ' - a'ap - pa 1 ' a), (10) 

1=1,2 1 ,=1,2 1 l 



where y^j and y,- are the pure dephasing rate and relaxation 
rate, respectively. We choose 7^,1 = 7^2 = J<p and j\ - yi - 7 
for identical DDMs. In simulation the coupling coefficient is 
chosen as g/2n =100 Mhz which has been experimentally re- 
alized f 27tl . To show the effect of pure dephasing clearly we do 
not assume the spin-echo technique is applied; therefore with 
the spin relaxation and dephasing times demonstrated above, 
the damping rates are translated to ~y<p/2n = 15.8 Mhz and 
yl2n = 0.16 Mhz. In simulation the photon number of the 
resonator is truncated at the value of 5. 



IV. QUANTUM COHERENCE VERSUS QUANTUM 
DISCORD 

Before present the results it is necessary to describe the 
quantification of quantum coherence and the optimization of 
the measurement to obtain the quantum discord. Through 
tracing over the states of the resonator, we get the re- 
duced density operator p^g = Tr f p, which describes the 



state of the composite system consisting of the two DDMs, 
and M, £ labels the individual DDM, respectively. To 
quantify the coherence of the composite system, we fo- 
cus on the the off-diagonal elements of p^<g in the basis 
{111)50, , |10>*b , lOD.^ , |00W. Yu and Song H have 
demonstrated the coherence in given basis can be quantified 
by the distance between the quantum state p^g and the nearest 
incoherent state. Thus the coherence of py\s can be measured 
by 

£W) = \\pms ~ tr*||i = Yj IP«4 (H) 

where || ■ ||i is the "entry wise" norm and cr* is the diagonal 
matrix with the identical diagonal elements of pw&. We note 
the measure T> includes the contribution of all off-diagonal 
elements of p^a, therefore it also contains the coherence of 
individual DDM J[ and S, which will have influence on the 
numerical results of the coherence dynamic for certain states. 
In the calculation of quantum discord, we choose the measure- 
ment set on & as the projectors {|^i) , (^Ih in which 
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FIG. 2. (Color online) Dynamics of quantum coherence and quantum 
correlations as a function of time. The solid (red), dashed (blue), 
and dotted (purple) curves correspond to the dynamics of discord, 
coherence, and EoF, respectively. The two DDM qubits are initially 
prepared in the state given in Eq. J 1 2b with (a) a 2 = 1/2, (b) a 2 = 
1/5, and (c) a 2 = 1/10. 



= cos6»|O) + e'' sin6)|l> and \if/ 2 ) = sin6> |0> - e'' cos 0|1>. 
The parameters 9 and <p are numerically varied from to 2n in 
order to seek the maximum of classical correlation C flU]. 

By numerical methods we obtain the dynamics of QC and 
QD in the two DDMs for different initial states. Firstly, we 
consider one kind of Bell-like state 



ar|01>*8+ Vl-a 2 |10>, 
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while the initial photon number in the TLR is assumed to 
be nonzero. For comparison we also calculate the entan- 
glement dynamics for the case a 2 = 1/2, using the EoF 
as the measurement of entanglement. For two qubits the 
EoF can be written in terms of Wootters' concurrence i34ll . 
E(t) = -/(f) log 2 f(t) - (1 - /(f)) log 2 (1 - /(f)), where f(t) = 

(1 + Vl - C(f) 2 )/2. The concurrence is expressed as C(t) = 
max{0, Ai - Aj— A3 — A4}, where A, (z = 1, 2, 3, 4) are the square 
roots of the eigenvalues of p^gcr,, ® o-yp^a-y ® cr y in decreas- 
ing order. The results are shown in Fig.|2fa)-|2jc). We observe 
that the QC, QD, and EoF all present oscillations and decays. 
For QC and QD there are two different types of amplitude 



of oscillation. The first type of amplitude undergoes decays 
after respective revivals, while the second type of amplitude 
increases and then decreases for a 2 - I and only decays for 
other values of a (an empty resonator also leads in this type of 
amplitude). For any pure state the QD and EoF coincide and 
for the maximally entangled state 2) = 1, thus for a 2 = 1/2 
the QC, QD, and EoF coincide initially. For non-maximally 
entanglement states the QC and QD certainly have different 
initial values. The EoF manifests the ESD and "sudden birth" 
phenomenons, which are the natural consequences under the 
combination action of the dissipation processes and the state 
transfer between DDMs and photons via the coupling to the 
TLR. Within short time intervals the EoF permanently van- 
ishes (see Fig. [2ja)) which is much earlier than the QC and 
QD, indicating the environment has greater impact on the EoF 
than QC and QD. Furthermore, the QD and QC vanish only at 
discrete instants and manifest the second type of amplitude of 
oscillation within the vanishing intervals of the EoF dynamics. 
This interesting feature shows some specific interaction with 
the environment does not induce the regular decay and revival 
of nonclassical correlations [8] (up to one vanishing instant), 
but leads to the revival of the nonclassical correlation as well 
as the coherence to a maximum then decays (two vanishing 
instants). In our model this feature is induced by the interac- 
tion with nuclear spins and low frequency fluctuations of the 
electrostatic bias, which are characterized by the pure dephas- 
ing term in the master equation in Eq. ( TTOb - The other feature 
is that the QC and QD exhibit homology to a considerable de- 
gree throughout for Bell-like states The QC and QD 
even vanish simultaneously. The discrepancy of the behaviors 
like the irregular shape of QD results from the dependance 
of QD on optimizing the measurement while QC is indepen- 
dent of measurement. We note that for Bell-like states pjyg 
has the X structure and the coherences of individual DDMs 
are always zero. Thus the QC for this case only captures the 
coherence of the joint two-qubit system. The results indicate 
that if there is no coherence in subsystems, the coherence dy- 
namic of the composite system shows the similar behavior of 
quantum discord, even with the absence of entanglement. 

Next we consider the initial states as the other kind of Bell- 
like state 
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while the TLR is assumed to be initially empty. The results 
are shown in Fig. |3ja)-[3jc). The behaviors of the QC, QD, 
and EoF are very different from those with the initial state 
\^i)jvB- The three dynamics present the similar behavior for 
f3 2 — 1 /2, but for f3 2 < 1/2 the oscillations and decays are 
quite complicated and exhibit respective trends. For some 
instants the QC and QD achieve a maximum while the EoF 
comes to a minimum, however despite the discrepancy of be- 
havior around the sudden changes of QD, the dynamics of QC 
and QD are similar, especially after a relative long time. The 
QD still manifests the robustness to the "sudden death" which 
is in agreement with previous works. Besides, unlike the first 
case the QC does not vanish within finite time intervals. 

At last we plot the dynamics of QC, QD, and EoF as a func- 
tion of time with the two DDMs prepared in a specific separa- 
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FIG. 3. (Color online) Dynamics of quantum coherence and quantum 
correlations as a function of time. The solid (red), dashed (blue), 
and dotted (purple) curves correspond to the dynamics of discord, 
coherence, and EoF, respectively. The two DDM qubits are initially 
prepared in the state given in Eq. J 1 3b with (a) B 2 = 1/2, (b)/? 2 = 1/5, 
and(c)y8 2 = 1/10. 
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FIG. 4. (Color online) Dynamics of quantum coherence and quantum 
correlations as a function of time. The solid (red), dashed (blue), 
and dotted (purple) curves correspond to the dynamics of discord, 
coherence, and EoF, respectively. The two DDM qubits are initially 
prepared in the separable state given in Eq. d 1 4b . 



ble state 

\<Pi)*s = 5(IO>y, + H>^)®(|0> s + |l> s ), (14) 

as shown in Fig.H Through the interaction with the photons 
in the resonator, the entanglement between the two DDMs 
is generated then oscillates until vanishing. Coupling to an 
empty resonator will lead to no entanglement but still the gen- 
eration of nonclassical correlations, in accordance with the 
results for the case of mediating by a common reservoir in 
Ref. One notices that the behaviors of QC and QD dif- 
fer widely, especially at the beginning, presenting substantial 
discrepancy between the initial separable state \<Pi)^ s and the 
Bell-like states; however after a relative long time, the be- 
haviors of QC and QD become similar. This is because for 
the initial state \ipi)^ B , there is no nonclassical correlations 
between the two DDMs but exist coherences in individual 
qubits, which are nonzero within short time intervals. In this 
case, the QC measure D captures not only the coherence of 
the composite system, but also the coherences in subsystems, 
while the QD measure <3 only quantifies the nonclassical cor- 
relations between subsystems. After a long time the cohereces 
in individual DDMs have been permanently lost due to the 
dissipations and the dynamics resembles the Bell-like cases, 
therefore the behaviors of QC and QD become similar after- 
wards. 

In experiments, the measurement of the entanglement be- 
tween distant DDMs is possible by the dispersive quantum 
nondemolition (QND) readout of qubits. In the dispersive 
regime, g/{a>d - u>o) <s 1, the QND readout of qubits can 
be realized by microwave irradiation of the resonator and then 
probing the transmitted or reflected photons. For example, 
by driving the resonator with a microwave pulse centered at 
the pulled frequency of cjq + g 1 l(tOd - u>o), the information 
about the state of the qubit is mostly stored in the number of 
transmitted photons. The scheme can be extended to the case 
of two or multiple qubits with different detunings (different 
DDMs) j26ll . and shows the feasibility of a joint measurement 
of entanglement in our system. 



V. SUMMARY AND PROSPECT 

In summary, we have numerically investigated the dynam- 
ics of quantum coherence and quantum correlations between 
two semiconductor DDMs separated by a distance and indi- 
rectly coupled via a transmission line resonator. The stronger 
impact on the EoF than QC and QD is verified, as well as the 
ESD phenomenon of EoF and the robustness of the QD to the 
sudden death, which implies the QD as a promising quantum 
information resource to process quantum computation in large 
scale solid systems. Furthermore, we find that the dephasing 
processes in our model induce a revival and decay of QC and 
QD within the vanishing interval of EoF for certain Bell-like 
states. We show that for initial Bell-like states, the QC and 
QD present similar behavior to a certain degree even with the 
absence of the entanglement, but for a initial specific separa- 
ble state the behaviors are totally different at first but become 
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similar after a long time. We find that the local coherences in 
individual qubits play an important role to the similarities and 
differences of the dynamics of QC and QD. The specific inter- 
nal relation between coherence and discord of bipartite mixed 
states may be revealed in an exactly solvable model while con- 
sidering a more proper coherence measure, ej* some measure 
connected with the localizable coherence 12911 . This inspires 
us to follow up in further research. 
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